REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1 21 5  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

06-11-1997  Paper 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

The  Influence  of  Combustion  Noise  on  Acoustic  Instabilities 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

V.S.  Burnely;  F.E.C.  Culick 

5d.  PROJECT  NUMBER 

3058 

5e.  TASK  NUMBER 

RF9A 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRS 

5  Pollux  Drive 

Edwards  AFBCA  93524-7048 

8.  PERFORMING  ORGANIZATION 
REPORT 

9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  (AFMC) 

AFRL/PRS 

5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 

10.  SPONSOR/MONITOR’S 
ACRONYM(S) 

11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-TP-1998-01 1 

12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


20020115  08/ 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 
PERSON 

Doug  Talley 

a.  REPORT 
Unclassified 

b.  ABSTRACT 
Unclassified 

c.  THIS  PAGE 
Unclassified 

A 

19b.  TELEPHONE  NUMBER 

(include  area  code) 

(661)  275-6174 

Standard  Form  298  (Rev.  8-98) 
Prescribed  by  ANSI  Std.  239.18 


The  Influence  of  Combustion  Noise  on  Acoustic  Instabilities 


V.  S.  Burnley* 

U.  S.  Air  Force  Research  Laboratory,  Edwards  Air  Force  Base,  California  93524 

and 

F.  E.  C.  Culickt 

California  Institute  of  Technology,  Pasadena,  California  91125 


Abstract 

Although  flows  in  combustors  contain  considerable  noise,  arising 
from  several  kinds  of  sources,  there  is  a  sound  basis  for  treating  or¬ 
ganized  oscillations  as  distinct  motions.  That  has  been  an  essential 
assumption  incorporated  in  virtually  all  treatments  of  combustion 
instabilities.  However,  certain  characteristics  of  the  organized  or  de¬ 
terministic  motions  seem  to  have  the  nature  of  stochastic  processes. 

For  example,  the  amplitudes  in  limit  cycles  always  exhibit  a  random 
character  and  even  the  occurrence  of  instabilities  seems  occasionally 
to  possess  some  statistical  features.  Analysis  of  nonlinear  coherent 
motions  in  the  presence  of  stochastic  sources  is  therefore  an  impor¬ 
tant  part  of  the  theory.  We  report  here  a  few  results  for  organized 
oscillations  in  the  presence  of  noise.  The  most  significant  deficiency 
of  this  work  is  that,  owing  to  the  low  level  of  current  understanding, 
the  stochastic  sources  of  noise  are  modeled  in  ad  hoc  fashion  and 
’Research  Engineer,  Propulsion  Directorate.  Member  AIAA. 

^Richard  L.  and  Dorothy  M.  Hayman  Professor  of  Mechanical  Engineering  and  Professor  of  Jet 
Propulsion.  Fellow  AIAA. 
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are  not  founded  on  a  solid  physical  basis  appropriate  to  combustion 
chambers. 
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Nomenclature 

=  speed  of  sound 

=  nonlinear  functional  arising  from  the  boundary  conditions 
=  forcing  function  of  the  nth  acoustic  mode 
=  nonlinear  functional  arising  from  the  conservation  equations 
=  wavenumber  of  the  nth  acoustic  mode 
=  unit  outward  normal  vector 
=  pressure 

=  magnitude  of  time-dependent  amplitude  rjn(t),  Eq.  (20) 

=  entropy 
=  velocity 

=  linear  growth  rate  of  the  nth  acoustic  mode 
=  time-dependent  amplitude  of  the  nth  acoustic  mode 
=  parametric  stochastic  excitation  associated  with  pressure 
=  parametric  stochastic  excitation  associated  with  velocity 
=  external  stochastic  excitation 
=  density 
=  noise  intensity 
=  correlation  time 

=  linear  frequency  shift  of  the  nth  acoustic  mode 
=  phase  of  the  time-dependent  amplitude  rjn(t),  Eq.  (20) 

=  mode-shape  of  the  nth  acoustic  mode 
=  frequency  of  the  nth  acoustic  mode 
=  vorticity 

=  acoustic  waves 
=  entropy  waves 
=  vorticity  waves 

=  nonlinear 
=  mean  quantity 
=  fluctuating  quantity 
=  time  derivative 
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1  Introduction 


Combustion  chambers  are  inherently  noisy  environments.  This  is  apparent  from  in¬ 
spection  of  the  power  spectra  of  pressure  records  from  test  firings,  as  well  as  from 
simply  listening  to  the  tests.  When  an  organized  oscillation,  i.e.,  a  combustion  in¬ 
stability,  is  present,  the  power  spectrum  exhibits  well-defined  peaks  in  addition  to 
background  noise  over  the  entire  range  of  frequencies.  Substantial  noise  sources  in 
rocket  motors  include  flow  separation,  turbulence,  and  combustion  processes.  It  is 
expected  that  the  presence  of  noise  will  affect  in  some  way  the  amplitudes  and  pos¬ 
sibly  the  qualitative  behavior  of  coherent  oscillations.  Nonetheless,  there  is  a  sound 
basis  for  treating  instabilities  as  distinct  motions,  and  much  progress  has  been  made 
using  this  assumption. 

There  are  two  general  types  of  coherent  oscillations  which  commonly  occur  in  com¬ 
bustion  chambers:  spontaneous  oscillations  and  pulsed  oscillations.  A  spontaneous 
oscillation  occurs  when  the  system  is  linearly  unstable.  As  a  result,  any  perturbation 
of  the  system  grows  exponentially  in  time.  Under  the  influence  of  nonlinear  effects, 
the  pressure  field  may  reach  a  periodic  motion,  or  limit  cycle.  This  type  of  instability 
is  also  known  as  an  intrinsic  instability  or  a  soft  excitation.  In  the  field  of  dynamical 
systems,  the  change  of  behavior  from  a  linearly  stable  steady  state  to  a  stable  path 
of  periodic  solutions  is  called  a  supercritical  bifurcation. 

A  pulsed  oscillation,  on  the  other  hand,  is  a  true  nonlinear  instability  of  a  lin¬ 
early  stable  system.  Small  perturbations  in  the  pressure  field  decay  exponentially 
to  zero,  while  larger  perturbations  may  lead  to  stable  or  unstable  periodic  motions. 
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Common  terminology  used  to  describe  this  type  of  oscillation  includes  triggered  in¬ 
stability,  hard  excitation,  and  subcritical  bifurcation.  Previous  works,  e.g.,  Yang  et 
al.,1  Paparizos  and  Culick,2  and  Jahnke  and  Culick,3  have  convincingly  shown  that 
nonlinear  gasdynamics  alone  does  not  contain  the  possibility  of  pulsed  oscillations. 
Part  of  the  purpose  of  the  present  investigation  is  to  determine  if  noise  processes 
could  be  responsible  for  the  occurrence  of  pulsed  oscillations. 

Culick  et  al.4  studied  the  influence  of  noise  on  combustion  instabilities,  but  only  for 
a  very  simple  case:  two  acoustic  modes  with  noise  present  only  in  the  first  mode.  In 
addition,  the  formulation  was  flawed,  and  the  form  of  the  resulting  noise  terms  is  not 
quite  correct.  Clavin  et  al.5  studied  the  influence  of  turbulence  on  instabilities  in  liquid 
rocket  motors.  Using  only  one  mode  in  the  analysis  and  third-order  nonlinearities,  it 
was  reported  that  the  inclusion  of  noise  can  lead  to  the  possibility  of  triggering.  It 
is  a  well-known  result  that  a  single  third-order  equation  may  produce  a  subcritical 
bifurcation.  When  more  acoustic  modes  are  considered,  this  may  not  be  the  case,  as 
demonstrated  by  Yang  et  al.1  for  third-order  gasdynamics.  Therefore,  the  results  of 
Clavin  et  al.  may  not  be  applicable  in  general. 

The  present  analysis  is  an  extension  of  the  previous  work  by  Culick  et  al.4  The 
approximate  analysis  used  here  has  been  covered  in  many  other  works,  so  only  a  brief 
overview  of  the  method  is  presented;  for  more  details,  see,  e.g.,  the  review  by  Culick.6 
Then,  we  will  decompose  the  flow  field  into  acoustic  and  non-acoustic  parts  based 
on  the  approach  followed  by  Chu  and  Kovasznay.7  That  tactic  produces  a  formal 
representation  of  the  noise  terms  in  the  governing  equations.  Finally,  we  will  simplify 
the  equations  in  order  to  study  more  conveniently  the  possible  influences  of  noise  on 
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combustion  instabilities. 


2  Development  of  the  Approximate  Analysis 

In  order  to  keep  the  analysis  as  general  as  possible,  the  formulation  begins  with  the 
conservation  equations  for  two-phase  flow.  These  equations  are  then  rewritten  in  an 
equivalent  form  for  a  single  medium  having  the  mass-averaged  properties  of  the  two 
phases.  Subsequently,  a  wave  equation  for  the  pressure  is  developed,  along  with  the 
corresponding  boundary  condition. 


vV  -  —  -  h 

V  P  a2  dt2  k 

n  ■  Vp'  =  -/ 


(1) 

(2) 


The  functions  h  and  /  are  linear  and  nonlinear  functions  of  the  pressure  and  velocity 
perturbations.  As  an  approximation,  these  perturbations  are  expanded  as  a  synthesis 
of  classical  acoustic  modes  with  time- varying  amplitudes, 

P'(r,t)  =pjr/Vn(t)ipn(r)  (3) 

71=1 

“>.*)  =  flW-VMr)  (4) 

n=l  lKn 

where  ipn  is  the  mode  shape  and  T]n(t)  is  the  time-dependent  amplitude  of  the  nth 
classical  acoustic  mode.  After  substituting  Eq.  (3)  in  Eq.  (1),  the  equations  are 
spatially  averaged,  resulting  in  a  system  of  ordinary  differential  equations  describing 
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the  amplitudes  of  the  acoustic  modes. 


+  “IVn  =  Fn  (5) 

where  ujn  =  akn  and 

__  2 

F»  =  -^{J'l’«hdV  +  fi>JdS]i  (6) 

Thus,  the  problem  is  reduced  to  solving  for  the  amplitudes,  r)n(t).  This  approach 
is  very  general  and  can  accommodate  all  damping  and  amplification  mechanisms. 
The  most  difficult  part  of  the  problem  is  in  the  identification  and  modeling  of  the 
important  physical  processes. 

3  Splitting  the  Unsteady  Flow  Field  Into  Acoustic,  Vortical,  and 
Entropic  Modes  of  Propagation 

Fluctuations  in  a  compressible  fluid  can  be  decomposed  into  three  types  of  waves: 
acoustic  waves,  vorticity  waves,  and  entropy  waves.  A  thorough  discussion  of  this 
idea  is  presented  by  Chu  and  Kovasznay.7  In  the  limit  of  small  amplitudes,  the  three 
types  of  waves  propagate  independently  in  a  uniform  mean  flow,  but  are  coupled 
when  the  mean  flow  is  non-uniform.8  For  example,  the  pressure  in  an  acoustic  wave 
is  changed  slightly  by  the  presence  of  a  vorticity  or  entropy  wave  if  the  mean  flow  is 
not  uniform.  Coupling  between  the  types  of  waves  may  also  occur  at  the  boundaries 
of  the  chamber. 
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Although  noise  is  detected  as  pressure  waves,  the  sources  of  noise  are  associated 
with  the  presence  of  vorticity  fluctuations  (e.g.,  turbulence,  flow  separation,  etc.)  and 
entropy  or  non-isentropic  temperature  fluctuations.  Therefore,  decomposing  the  un¬ 
steady  flow  field  into  the  three  types  of  waves  allows  both  noise  and  acoustic  insta¬ 
bilities  to  be  handled  in  the  same  analytical  framework  discussed  in  Sec.  1.  The 
contributions  from  vorticity  and  entropy  waves  appear  as  additional  force  terms  on 
the  right-hand  side  of  the  acoustic  equation. 

Following  the  analysis  of  Chu  and  Kovasznay,7  the  thermodynamic  and  kinematic 
variables  can  be  written  as  sums  of  fluctuations  in  the  three  waves  as  follows. 


P' =p'a+Pa+P's  (7) 

ft'  =  fya  +  Q'n  +  n;  (8) 

s'  =  s'a  +  Sfi  +  s'  (9) 

u'  =  u'a  +  u'u  +  u's  (10) 


In  general,  all  of  the  fluctuations  will  be  nonzero,  but  not  all  terms  are  of  the  same 
order.  If  we  restrict  the  analysis  to  small  amplitude  motions,  the  three  waves  have 
the  following  characteristics:7 

(  )a  acoustic  waves:  pressure  and  velocity  fluctuations,  no  entropy  change 
(  )n  vorticity  waves:  velocity  fluctuations,  no  pressure  or  entropy  changes 
(  )<,  entropy  waves:  entropy  and  velocity  fluctuations,  no  pressure  change 


7 


Thus,  to  zeroth-order,  the  fluctuations  in  the  three  waves  are  given  by- 


P'  =Pa 

(11) 

^  = 

(12) 

s'  =  s'3 

(13) 

u'  =  u'a  +  u'n  +  u's 

(14) 

An  equation  for  the  density  fluctuation  is  obtained  by  expanding  the  formula  for  the 
entropy  of  a  perfect  gas. 


P  7  P 


(15) 


As  an  approximation  to  the  acoustic  pressure  and  velocity  perturbations,  we  will 
once  again  use  a  superposition  of  the  classical  acoustic  modes  so  that 


0°  oo  ■  (f\ 

Pa  =  P  E  Vn{t)Mr)  u'a  =  £  ' ^Mr) 

n= 1  n=l  7 Kn 

Substitution  of  Eq.  (11)  in  the  left-hand  side  of  the  nonlinear  wave  equation  (1), 
followed  by  application  of  Galerkin’s  method,  leads  to  a  set  of  coupled  nonlinear 
oscillator  equations. 


V n  +  ulvn  =  Fn 


(16) 
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where8 

-  ^f-Fn  =  ph  +  +  ph  +  z^h 

a 1  az  cr 

r  f)uf  r  1  dVf 

+l^-hds~f <i7) 

and 

h  =  J  (u  •  Vu'  +  u'  ■  Vu)  •  VipndV  J2  =  J (72/V  •  u  +  uV  •  p')ij)ndV 

i3  =  J  (u'  •  v«i'  +  j  •  v^n<iy  i,  =  ^  J  (7?'v  •  u’  +  u'  • 

In  the  original  development  of  the  approximate  analysis,  the  zeroth-order  approxi¬ 
mations  for  the  pressure  and  velocity  were  used  to  evaluate  Fn.  The  same  idea  will 
be  applied  here,  although  additional  contributions  to  the  velocity  fluctuation  from 
coupling  to  vorticity  and  entropy  waves  will  be  included.  Once  these  quantities  are 
substituted  in  the  right-hand  side,  the  set  of  forced  oscillator  equations  eventually 
takes  the  general  form  (see  Burnley9  for  details) : 

OO  OO 

fjn  “I"  ^  ^  ^  ^  [^nijViVj  "1“  -^nijViVj] 

2=1  j  —  1 

OO 

+  (-^n) other  +  H  ICiiVi  +  ZniVi]  +  -n  (18) 

i=l 

This  system  of  equations  is  very  complex  with  many  free  parameters.  For  in¬ 
stance,  if  we  truncate  the  system  to  N  modes,  there  are  2  N  linear  parameters  and  an 
additional  2N2  +  N  unknown  functions.  In  order  to  simplify  the  equations  somewhat, 
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we  will  therefore  neglect  cross-coupling  terms  in  and  i.e.,  terms  with  n  ±  i. 
These  terms  may  turn  out  to  be  important,  but  neglecting  them  will  allow  for  easier 
initial  computation  of  results  and  will  suffice  for  the  purposes  here.  The  simplified 
set  of  equations  is 


00  00 

Vn  “t"  ^nVn  ~  20!n f]n  +  2uJn6nT}n  ^  ^  T  Bnij Tji Tjj ] 

i=  1  j=l 

+  (^n)  other  +  ZnitfVn  +  €n(t)Vn  +  Hn(t)  (19) 

4  Modeling  of  the  Stochastic  Sources 

The  problem  has  now  been  reduced  to  solving  Eq.  (19)  for  the  time-dependent  am¬ 
plitudes  r]n(t)  which  are  subjected  to  additive  and  multiplicative  noise.  The  source 
terms  £”(i),  £n(t),  and  3n(t)  represent  stochastic  processes  of  some  sort  and  are  re¬ 
sponsible,  in  this  formulation,  for  the  background  noise  found  in  the  power  spectra 
of  test  firings.  The  problem  of  modeling  these  processes,  however,  remains.  This 
requires  specification  of  both  the  spatial  and  temporal  distribution  of  the  velocity 
and  the  entropy.  At  the  present  time,  no  models  exist  for  these  fluctuations. 

There  are  several  other  paths  that  can  be  followed  at  this  point  which  include 
obtaining  approximate  representations  for  the  velocity  and  entropy  fluctuations  based 
on  experimental  data  or  numerical  simulations.  The  approach  that  will  be  taken 
here  is  to  assume  forms  for  the  source  terms  which  are  based  on  observations  of 
experiments.  By  inspection  of  the  pressure  traces  of  test  firings,  it  is  apparent  that 
the  stochastic  processes  in  real  systems  are  broadband  with  very  small  correlation 
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times*,  tc.  The  limit  rc  -*  0  represents  a  delta  correlated  process,  i.e.,  a  process 
which  is  totally  uncorrelated  with  itself.  It  is  thus  interesting  to  study  this  limiting 
case  and  assume  that  the  stochastic  terms  are  represented  by  white  noise. 

The  definition  of  a  white  noise  process  is  a  process  whose  spectral  density  is  flat, 
i.e.,  all  frequencies  are  present  at  the  same  amplitude.  Although  such  a  process  cannot 
occur  in  a  real  system,  white  noise  can  be  a  very  useful  tool  for  studying  real  processes 
which  have  very  small  correlation  times  compared  to  the  macroscopic  times  of  the 
system.  This  is  true  of  the  random  processes  and  systems  of  interest.4  Therefore,  we 
will  approximate  £,  £„,  and  5„  by  mutually  independent  white  noise  processes  with 
zero  mean  values  and  intensities  denoted  by  <r^  ,  and  cn . 

An  example  of  a  simulation  with  white  noise  excitations  is  presented  in  Fig.  1. 
A  sample  trace  of  the  pressure  at  the  head  end  of  the  combustion  chamber  is  shown 
along  with  the  corresponding  normalized  spectrum.  Inspection  of  the  spectrum  shows 
peaks  at  distinct  frequencies  which  are  associated  with  an  acoustic  instability,  along 
with  broadband  background  noise.  This  is  characteristic  of  actual  test  data  of  a  case 
when  an  instability  is  present.  For  the  simulation  shown  in  Fig.  1,  as  well  as  all 
other  results  presented  here,  the  values  given  in  Table  1  will  be  used  for  the  linear 
parameters  of  the  system. 


*The  correlation  time  is  the  time  above  which  the  autocorrelation  function  is  zero.  This  is  a 
measure  of  the  dependence  of  the  process  on  its  past. 


11 


n 

1 

2 

3 

4 

an(sec~l) 

free 

-324.8 

-583.6 

-889.4 

8n  (rad/sec) 

12.9 

46.8 

-29.3 

-131.0 

Table  1:  Linear  growth  rates  and  frequency  shifts 


5  Results 

Since  we  are  interested  in  nondeterministic  systems,  it  is  natural  to  use  the  probability 
density  functions  of  the  amplitudes  of  acoustic  modes  to  investigate  the  dynamics  of 
the  system.  In  particular,  we  will  use  the  amplitude  rn  which  is  defined  in  the  following 
equation. 

Vn(t)  =  rn(t)  cos(unt  +  <j>n{t))  (20) 

Using  rn  instead  of  r]n,  any  quantitative  changes  which  occur  in  the  solution  will  be 
more  apparent  since  rn  is  a  measure  of  the  magnitude  of  the  oscillation.  Note  that 
this  change  of  variables  from  ( rjn ,  f]n)  to  ( rn ,  (fin)  will  cause  a  shift  in  the  expected 
value.  This  is  obvious  from  the  definition  of  rn,  which  is  defined  as  yrft  -f  rj2  Jul 
A  Monte-Carlo  method  will  be  used  to  obtain  an  approximation  to  the  probability 
density  functions  of  the  acoustic  amplitudes.  In  this  method,  a  series  of  numerical  “ex¬ 
periments”  is  conducted,  usually  in  the  same  manner  that  one  would  conduct  actual 
experiments.  After  the  flow  field  has  become  well-developed  (say  after  1000  periods 
of  the  fundamental  mode  or  so),  the  amplitudes  of  the  acoustic  modes  are  sampled. 
The  results  are  then  used  to  construct  histograms  which,  after  normalization,  approx- 
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imate  the  instantaneous  probability  density  functions  of  the  modal  amplitudes.  The 
approximation  becomes  better  as  the  number  of  experiments  is  increased. 

In  the  current  study,  each  Monte-Carlo  simulation  will  consist  of  10000  numerical 
experiments.  The  linear  parameters  will  be  fixed  throughout  a  series  of  experiments, 
while  the  initial  conditions  for  the  simulations  will  be  varied  systematically.  In  par¬ 
ticular,  a  square  initial  pulse  which  is  nonzero  from  0  <  x/L  <  .25  will  be  used.  The 
size  of  the  pressure  pulse  p'/p  will  be  varied  from  0  to  .2  linearly  so  as  to  include  most 
likely  values. 

For  the  initial  results,  we  will  use  the  simplest  possible  set  of  equations.  The 
system  will  be  truncated  to  two  modes  with  noise  included  explicitly  only  in  the 
fundamental  mode.  We  will  also  assume  initially  that  the  only  nonlinear  contributions 
are  associated  with  gasdynamics.  These  simplifications  were  used  by  Culick  et  al.4 
and  will  allow  the  effects  of  each  type  of  stochastic  excitation  to  be  determined.  We 
will  relax  these  simplifications  later  in  this  section  when  nonlinear  contributions  from 
combustion  are  also  included. 

5.1  The  Effects  of  an  External  Excitation 

The  term  Ei  is  an  external  excitation,  i.e.,  it  does  not  depend  on  the  current  state  of 
the  system.  As  a  result,  this  type  of  excitation  does  not  change  the  qualitative  behav¬ 
ior  of  the  system  from  that  of  the  original  deterministic  system,10  i.e.,  no  matter  how 
large  the  intensity,  the  addition  of  an  external  excitation  does  not  affect  the  locations 
or  types  of  attractors  in  the  system.  Instead,  an  external  excitation  provides  only  a 
random  perturbation  from  these  states;  the  dynamics  of  the  system  are  constantly 
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acting  to  bring  a  stable  system  back  to  a  state  of  equilibrium. 

To  understand  better  the  effects  of  an  external  excitation,  it  is  useful  to  treat  the 
deterministic  case  as  the  limit  of  of  0.  For  a  deterministic  system  with  nonlinear 
contributions  from  gasdynamics  only,  the  stationary  probability  density  functions  will 
be  delta  functions  in  terms  of  rn.  For  example,  the  probability  density  functions  of  a 
linearly  stable  system  without  random  perturbations  will  be  delta  functions  at  rn  =  0 
for  all  n.  When  an  external  excitation  is  introduced,  the  effect  is  to  shift  the  mean 
value  of  the  amplitude  rn  to  a  nonzero  value  (due  to  the  change  of  variables)  and 
increase  the  variance  of  the  oscillation  such  that  a  broader  range  of  values  is  likely. 
This  broadening  effect  can  be  seen  in  Figs.  2  and  3  for  linearly  stable  and  linearly 
unstable  systems,  respectively. 

To  demonstrate  the  effect  of  noise  intensity  on  the  probability  density,  the  value 
used  in  Fig.  3  was  doubled,  and  the  result  is  plotted  in  Fig.  4.  As  the  intensity  of  the 
noise  increases,  the  variance  of  the  probability  density  function  increases  so  that  a 
larger  range  of  modal  amplitudes  is  likely.  The  mean  value  of  the  amplitude,  however, 
remains  unchanged. 

Another  interesting  aspect  of  an  external  excitation  can  be  seen  if  we  look  at 
a  sample  pressure  history  and  spectral  density  for  a  linearly  stable  system  with  an 
external  excitation.  At  first  inspection  of  Fig.  5,  it  might  appear  that  the  system  is 
linearly  unstable  due  to  the  very  low  amplitude  oscillation.  The  data  shown  in  Fig.  5 
is,  however,  for  a  stable  system.  It  happens  that  an  external  excitation,  though  ran¬ 
dom,  will  excite  the  natural  acoustic  modes  of  the  chamber  such  that  a  low  amplitude 
fluctuation  is  present.  The  energy  provided  by  the  external  excitation  Ex  is  most  eas- 
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ily  transferred  into  the  fundamental  acoustic  mode  because  this  excitation  appears 
in  the  equation  for  the  fundamental  mode.  Therefore,  the  system  oscillates  at  the 
fundamental  frequency  of  the  chamber.  However,  in  terms  of  dynamical  systems  the¬ 
ory,  the  system  is  still  linearly  stable,  and  the  dynamics  of  the  system  always  act  to 
bring  the  chamber  back  to  a  state  of  rest.  Recently,  in  analyzing  data  similar  to  this, 
Malhotra  and  Flandro11  suggested  that  the  system  was  linearly  unstable  due  to  the 
presence  of  a  very  low  amplitude  oscillation.  We  believe  that  the  system  being  ana¬ 
lyzed  by  Flandro  was  actually  linearly  stable  and  that  the  low  amplitude  oscillation 
was  actually  caused  by  the  noise  processes  which  are  always  present  in  real  systems. 

5.2  The  Effects  of  a  Noisy  Linear  Growth  Rate 

Unlike  the  external  excitations  covered  in  the  previous  section,  the  parametric  ex¬ 
citations  ?)i  and  ^rji  depend  upon  the  current  state  of  the  system.  In  a  more 
illuminating  form,  the  system  (19)  can  be  rewritten  as 

Vi  +UJirh  =  2  ^ai  +  771  +  2ui  (di  +  — rji  +  (Fi)nl  +  Hi (t) 

V2  +  ojIt]2  =  2a2T]2  +  2uj202r)2  +  (F2)NL  (21) 

By  inspection  of  the  above  system,  it  is  easy  to  see  that  £[  is  a  random  perturbation 
of  the  linear  growth  rate  of  the  first  acoustic  mode.  Similarly,  £1  is  a  random  per¬ 
turbation  of  the  linear  frequency  shift.  In  this  section,  we  will  study  the  effects  of  a 
noisy  linear  growth  rate  on  the  dynamics  of  the  system,  while  the  effects  of  a  noisy 
linear  frequency  shift  will  be  handled  in  the  next  section. 
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There  are  two  main  effects  caused  by  the  noisy  linear  growth  rate  ft.  The  first  is 
a  result  of  approximating  a  real  noise  process  by  white  noise.  Since  no  real  process  is 
truly  white,  the  system  will  have  a  small  but  finite  memory,  i.e.,  there  will  be  some 
correlation  between  the  noise  and  the  system.  This  correlation  is  taken  into  account 
by  the  Stratonovich  representation  through  its  definition  of  the  stochastic  integral.10 
As  a  result,  the  linear  growth  rate  is  increased  to  an  apparent  value  given  by 

/  \  2 

(^l) apparent  =  Oi\  +  f  J  (22) 

Thus,  one  effect  of  a  noisy  linear  growth  rate  is  to  shift  the  bifurcation  diagram  by 
an  amount  proportional  to  the  square  of  the  intensity  of  the  noise.  This  is  known  as 
noise-induced  drift.10  The  second  effect  of  the  noisy  linear  growth  rate  is  similar  to 
the  effect  of  an  external  excitation.  It  is  basically  a  disorganizing  effect  which  tends 
to  spread  the  peak  of  the  probability  density  functions  about  the  mean  value. 

This  type  of  parametric  excitation  was  studied  by  Horsthemke  and  Lefever10  on 
a  first-order  nonlinear  equation.  In  that  study,  the  Verhulst  model  with  a  noisy 
growth  rate  was  shown  to  have  two  transition  points  at  which  the  probability  density 
function  changes  qualitatively.  Thus,  three  ranges  of  linear  growth  rate  which  produce 
three  different  types  of  probability  density  function  were  found.  Three  qualitatively 
different  types  of  probability  density  functions  are  also  found  in  our  system.  However, 
it  is  difficult  to  predict  the  exact  location  of  the  transitions  due  to  the  highly  nonlinear 
nature  of  the  system  of  equations. 

For  highly  stable  systems,  i.e.,  ai  ^  0  sec  1 ,  the  attractor  at  rn  =  0  is  so  strong 
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that  the  solution  will  always  be  drawn  back  to  this  steady  state.  Once  this  state 
is  reached,  the  parametric  excitation  no  longer  affects  the  system.  As  a  result,  the 
stationary  probability  density  functions  for  highly  stable  systems  will  be  delta  func¬ 
tions  at  rn  =  0  for  all  n.  For  mildly  stable  and/or  unstable  systems,  a  second  type  of 
probability  density  function  occurs.  While  the  most  probable  value  is  still  zero,  the 
mean  value  is  not.  This  type  of  probability  density  function  was  shown  previously  in 
the  study  by  Culick  et  al.4  An  example  is  shown  here  in  Fig.  6  for  ax  =  -8  sec-1. 

The  final  type  of  probability  density  occurs  in  systems  which  are  highly  unstable. 
In  that  case,  the  linear  growth  rate  ax  is  so  large  that  it  is  unlikely  that  the  noise 
will  be  strong  enough  for  a  sufficient  amount  of  time  to  drive  the  solution  back  to 
the  trivial  state.  Therefore,  both  the  mean  value  and  the  most  probable  value  are 
nonzero.  Figure  7  shows  an  example  of  such  a  probability  density  function  for  ax  =  25 
sec"1.  Note  that  the  mean  value  is  shifted  slightly  from  the  deterministic  value  of 
0.07.  This  is  an  example  of  noise-induced  drift. 

For  the  deterministic  system  with  nonlinear  contributions  from  gasdynamics  only, 
two  qualitatively  different  regions  are  found:  one  region  of  stable  steady  states  and 
one  region  of  stable  limit  cycles.  When  noise  of  a  parametric  nature  is  added,  three 
distinct  regions  arise  because  rn  =  0  is  a  stationary  point  of  the  system  for  all  values  of 
and  CT^-10  However,  if  an  external  excitation  is  also  included  in  the  system,  rn  =  0 
is  no  longer  a  stationary  state  (although  it  is  still  an  attractor),  and  only  two  different 
regions  are  found. 
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5.3  The  Effects  of  a  Noisy  Linear  Frequency  Shift 

The  effects  of  a  noisy  linear  frequency  shift  are  very  similar  to  the  effects  of£\  so  only 
a  brief  discussion  is  necessary.  Since  is  also  a  parametric  excitation,  a  stationary 
state  occurs  once  again  for  rn  =  0,  and  three  regions  of  distinct  types  of  probability 
density  functions  are  produced.  The  three  types  are  qualitatively  similar  to  those  in 
the  previous  section,  although  the  transition  points  may  occur  at  different  values  of 
^  sample  probability  density  function  for  a  linearly  unstable  system  with  a  noisy 
linear  frequency  shift  is  plotted  in  Fig.  8. 

5.4  The  Effects  of  Noise  and  Nonlinear  Combustion 

In  the  previous  sections  involving  nonlinear  contributions  from  gasdynamics  only,  no 
cases  were  found  which  are  consistent  with  the  qualitatively  different  behaviour  of 
triggering,  i.e.,  nonlinear  instabilities  in  a  linearly  stable  system.  As  no  examples  of 
triggering  were  found  for  the  case  of  noise  and  nonlinear  gasdynamics,  a  model  of 
nonlinear  combustion  which  has  previously  been  shown12  to  provide  the  possibility  of 
triggering  will  be  included  in  the  analysis  in  order  to  determine  the  effects  of  noise  on 
a  deterministic  system  which  has  multiple  stable  stationary  states  for  the  same  value 

Of  CH\. 

The  threshold  velocity  model  is  an  ad  hoc  model  which  is  based  on  the  idea  of 
velocity  coupling  with  a  threshold  value  below  which  the  effects  of  nonlinear  combus¬ 
tion  are  not  felt.  Threshold  effects  have  been  observed  in  experimental  investigations 
of  velocity  coupling.  Ma  et  al.13  found  a  threshold  acoustic  velocity  above  which  the 
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mean  mass  flux  increased  linearly  with  the  Reynold’s  number  of  the  acoustic  fluctu¬ 
ations.  Below  the  threshold  value,  the  mean  mass  flow  was  approximately  constant. 
It  was  determined  that  the  increased  mass  flux  is  a  result  of  increased  heat  transfer 
to  the  surface  after  transition  to  turbulent  flow  had  occurred. 

The  form  of  the  threshold  velocity  model  is  given  by 

fa'  =  rhRvcF(u')  (23) 

where  the  function  F(u')  is  shown  in  Fig.  9.  This  function  introduces  a  dead  zone  in 
which  the  nonlinear  contributions  from  combustion  do  not  affect  the  system.  When 
the  amplitudes  of  oscillations  become  larger  than  the  chosen  threshold  value  ut,  the 
nonlinear  effects  are  then  felt.  The  threshold  velocity  model  can  be  included  in  the 
analysis  using  the  term  (F *)other>‘  the  details  of  this  are  given  in  Burnley.9 

When  nonlinear  combustion  is  added  to  the  stochastic  system  (19),  the  resulting 
probability  density  functions  can  be  quite  different,  as  one  might  expect.  In  a  previous 
investigation,  we  have  shown  that  the  threshold  velocity  model  can  produce  regions  of 
possible  triggering  in  which  two  stable  solutions  exist  simultaneously.  In  a  stochastic 
system,  this  corresponds  to  a  bimodal  probability  density  function  such  that  there  is 

a  high  probability  of  low  and  high  amplitudes  and  a  low  probability  of  intermediate 
values. 

For  the  parametric  values  in  the  threshold  velocity  model,  we  will  use  a  nondi- 
mensional  threshold  velocity  of  ut/a  =  0.03  and  Rvc  =  7.8.  Using  these  values,  the 
bifurcation  diagram  for  the  deterministic  system  is  shown  in  Fig.  10.  This  diagram 
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will  be  useful  in  the  discussion  of  results.  In  addition,  the  following  values  were  chosen 
for  the  intensities  of  the  stochastic  sources:  =  0.005  sec-3/2,  cr£  =  0.025  sec-1/2, 

and  cxn  —  0.0005  sec  3//2,  for  n  =  1, 2.  The  parameter  oc\  will  be  varied  while  all  other 
parametric  values  remain  fixed.  By  changing  this  parameter,  we  will  demonstrate  a 
variety  of  the  possible  forms  of  the  probability  density  functions. 

From  inspection  of  Fig.  10,  we  see  that  the  region  of  possible  triggering  begins  at 
approximately  a\  =  —30  sec  1  for  the  deterministic  system.  Below  this  value,  the 
deterministic  system  is  stable  to  any  size  perturbation.  To  illustrate  the  effect  of  noise 
on  such  a  system,  a  linear  growth  rate  of  -35  sec-1  was  chosen.  Figure  12a  shows 
the  resulting  probability  density  function  for  the  first  acoustic  mode.  For  this  case, 
the  attractor  of  the  deterministic  system,  i.e.,  the  trivial  steady  state,  is  so  strong 
that  the  amplitudes  never  reach  large  values.  Therefore,  the  parametric  excitations, 
and  £7 1,  have  a  very  small  effect  on  the  system.  Most  of  the  noise  contribution 
is  a  result  of  the  external  excitations,  Sn. 

As  the  linear  growth  rate  is  increased  to  a  value  above  -30  sec-1,  we  enter  the  re¬ 
gion  of  possible  triggering  for  the  deterministic  system  where  an  additional  attractive 
state  is  present.  Three  values  of  ft]  were  chosen  in  order  to  show  how  the  probability 
density  of  the  fundamental  mode  changes  throughout  this  region.  As  the  value  of  ax 
is  varied,  the  regions  of  attraction  of  the  stationary  states  will  change.  This  will  have 
a  noticeable  effect  on  the  probability  density  functions. 

The  first  value  of  ft]  is  very  close  to  the  lower  boundary  of  the  region  of  possible 
triggering.  For  a  value  of  ft]  =  —25  sec  1 ,  Fig.  12b  shows  the  probability  density 
function  of  the  fundamental  mode.  The  low  amplitude  attractive  state  is  dominant 
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Figure  1:  Sample  pressure  trace  and  spectrum  for  a  simulation 


Figure  2:  The  influence  of  an  external  excitation  of  the  fundamental  mode  only  on  a 
linearly  stable  system;  2  modes,  af  =  0.0005  sec-3/2,  =  —25  sec-1 
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Figure  4:  The  influence  of  increasing  the  intensity  of  an  external  excitation  of  the 
fundamental  mode  only;  2  modes,  of  =  0.001  sec-3/2,  =  25  sec-1 
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Figure  5:  Sample  pressure  trace  and  spectrum  for  a  simulation  of  a  stable  system 
with  an  external  excitation;  2modes,  erf  =  0.0005  sec-3/2,  oi\  =  —25  sec-1 


Figure  6:  Probability  density  function  for  a  system  with  a  noisy  linear  growth  rate  in 
the  fundamental  mode  only;  2  modes,  erf  =  0.01  sec-3/2,  au  =  —8  sec-1 
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Figure  7:  Probability  density  function  for  a  system  with  a  noisy  linear  growth  rate  in 
the  fundamental  mode  only;  2  modes,  of  =  0.005  sec-3/2,  a\  =  25  sec-1 


Figure  8:  Probability  density  function  for  a  system  with  a  noisy  linear  frequency  shift 
in  the  fundamental  mode  only;  2  modes,  of  =  0.025  sec-1/2,  oi\  =  25  sec-1 
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Figure  9:  Threshold  velocity  function;  ut/a  =  0.03 


Figure  10:  Bifurcation  diagram  for  the  deterministic  system;  threshold  velocity  model, 
ut/a  =  0.03,  RyC  =  7.8,  4  modes 


List  of  Figures 


4 


* 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


Sample  pressure  trace  and  spectrum  for  a  simulation . 

The  influence  of  an  external  excitation  of  the  fundamental  mode  only 
on  a  linearly  stable  system;  2  modes,  erf  =  0.0005  sec-3/2,  a.\  =  —25 

sec-1  . 

The  influence  of  an  external  excitation  of  the  fundamental  mode  only 
on  a  linearly  unstable  system;  2  modes,  erf  =  0.0005  sec-3/2,  ot\  —  25 

sec-1  . 

The  influence  of  increasing  the  intensity  of  an  external  excitation  of 
the  fundamental  mode  only;  2  modes,  erf  =  0.001  sec-3/2,  =  25  sec' 

Sample  pressure  trace  and  spectrum  for  a  simulation  of  a  stable  system 
with  an  external  excitation;  2modes,  erf  =  0.0005  sec-3/2,  a.\  —  —25 

sec-1  . 

Probability  density  function  for  a  system  with  a  noisy  linear  growth 
rate  in  the  fundamental  mode  only;  2  modes,  of  —  0.01  sec-3/2,  a\  = 

—8  sec-1 . 

Probability  density  function  for  a  system  with  a  noisy  linear  growth 
rate  in  the  fundamental  mode  only;  2  modes,  of  =  0.005  sec-3/2, 

cci  =  25  sec-1 . 

Probability  density  function  for  a  system  with  a  noisy  linear  frequency 
shift  in  the  fundamental  mode  only;  2  modes,  of  =  0.025  sec-1/2, 

d\  =  25  sec-1 . 

Threshold  velocity  function;  ut/a  —  0.03 . 

Bifurcation  diagram  for  the  deterministic  system;  threshold  velocity 

model,  ut/a  =  0.03,  Rvc  =  7.8,  4  modes . 

Bifurcation  diagram  for  the  deterministic  system;  2  modes . 

Probability  density  function  of  the  first  acoustic  mode  for  various  values 
of  au;  threshold  velocity  model,  ut/a  =  0.03,  Rvc  =  7.8,  4  modes  .  .  . 
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